Next Generation Sequencing, and Development of a Pipeline as a Tool for the Detection and Discovery of Citrus Pathogens to Facilitate Safer Germplasm Exchange.

Citrus is affected by many diseases, and hence, the movement of citrus propagative materials is highly regulated in the USA. Currently used regulatory pathogen detection methods include biological and laboratory-based technologies, which are time-consuming, expensive, and have many limitations. There is an urgent need to develop alternate, rapid, economical, and reliable testing methods for safe germplasm exchange. Citrus huanglongbing (HLB) has devastated citrus industries leading to an increased need for germplasm exchanges between citrus growing regions for evaluating many potentially valuable hybrids for both HLB resistance and multilocational performance. In the present study, Next-Generation Sequencing (NGS) methods were used to sequence the transcriptomes of 21 test samples, including 15 well-characterized pathogen-positive plants. A workflow was designed in the CLC Genomics Workbench software, v 21.0.5 for bioinformatics analysis of the sequence data for the detection of pathogens. NGS was rapid and found to be a valuable technique for the detection of viral and bacterial pathogens, and for the discovery of new citrus viruses, complementary to the existing array of biological and laboratory assays. Using NGS methods, we detected beet western yellows virus, a newly reported citrus virus, and a variant of the citrus yellow vein-associated virus associated with the "fatal yellows" disease.

Variation in microbial feature perception in the Rutaceae family with immune receptor conservation in citrus.

Although much is known about the responses of model plants to microbial features, we still lack an understanding of the extent of variation in immune perception across members of a plant family. In this work, we analyzed immune responses in Citrus and wild relatives, surveying 86 Rutaceae genotypes with differing leaf morphologies and disease resistances. We found that responses to microbial features vary both within and between members. Species in 2 subtribes, the Balsamocitrinae and Clauseninae, can recognize flagellin (flg22), cold shock protein (csp22), and chitin, including 1 feature from Candidatus Liberibacter species (csp22CLas), the bacterium associated with Huanglongbing. We investigated differences at the receptor level for the flagellin receptor FLAGELLIN SENSING 2 (FLS2) and the chitin receptor LYSIN MOTIF RECEPTOR KINASE 5 (LYK5) in citrus genotypes. We characterized 2 genetically linked FLS2 homologs from "Frost Lisbon" lemon (Citrus ×limon, responsive) and "Washington navel" orange (Citrus ×aurantium, nonresponsive). Surprisingly, FLS2 homologs from responsive and nonresponsive genotypes were expressed in Citrus and functional when transferred to a heterologous system. "Washington navel" orange weakly responded to chitin, whereas "Tango" mandarin (C. ×aurantium) exhibited a robust response. LYK5 alleles were identical or nearly identical between the 2 genotypes and complemented the Arabidopsis (Arabidopsis thaliana) lyk4/lyk5-2 mutant with respect to chitin perception. Collectively, our data indicate that differences in chitin and flg22 perception in these citrus genotypes are not the results of sequence polymorphisms at the receptor level. These findings shed light on the diversity of perception of microbial features and highlight genotypes capable of recognizing polymorphic pathogen features.

The Development and Evaluation of Insect Traps for the Asian Citrus Psyllid, Diaphorina citri (Hemiptera: Psyllidae), Vector of Citrus Huanglongbing.

Citrus huanglongbing (HLB) is a severe problem for citrus cultivation. The disease management programs benefit from improved field tools suitable for surveying the ACP vector (Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Psyllidae)) and the associated pathogen. In the present study, we utilize three-dimensional (3D) printers and design tools to develop traps that can capture and preserve ACPs. Three novel, 3D-printed traps were designed and evaluated: stem trap, and cylinder traps 1 and 2. The traps and yellow sticky cards were deployed weekly for 8 months in 2 non-commercial citrus groves in Florida; in California, the traps were evaluated for 12 months in field cages and 4 citrus groves. The stem traps captured lower numbers of ACPs at all experimental sites compared to the cylinder traps. Capture rates in the cylinder traps were comparable to the sticky trap, making the device a viable tool for monitoring field ACPs. The two main advantages of using the reusable 3D traps over standard methods of ACP and HLB surveys include dynamic sampling that can be conducted year-round and the capture of ACPs that can be preserved and tested. Improved trapping may facilitate quick management decisions and mitigate HLB.

An Improved Reference Gene for Detection of "Candidatus Liberibacter asiaticus" Associated with Citrus Huanglongbing by qPCR and Digital Droplet PCR Assays.

Citrus huanglongbing (HLB) disease associated with the 'Candidatus Liberibacter asiaticus' (CLas) bacterium has caused significant financial damage to many citrus industries. Large-scale pathogen surveys are routinely conducted in California to detect CLas early in the disease cycle by lab-based qPCR assays. We have developed an improved reference gene for the sensitive detection of CLas from plants in diagnostic duplex qPCR and analytical digital droplet PCR (ddPCR) assays. The mitochondrial cytochrome oxidase gene (COX), widely used as a reference, is not ideal because its high copy number can inhibit amplification of small quantities of target genes. In ddPCRs, oversaturation of droplets complicates data normalization and quantification. The variable copy numbers of COX gene in metabolically active young tissue, greenhouse plants, and citrus relatives suggest the need for a non-variable, nuclear, low copy, universal reference gene for analysis of HLB hosts. The single-copy nuclear gene, malate dehydrogenase (MDH), developed here as a reference gene, is amenable to data normalization, suitable for duplex qPCR and ddPCR assays. The sequence of MDH fragment selected is conserved in most HLB hosts in the taxonomic group Aurantioideae. This study emphasizes the need to develop standard guidelines for reference genes in DNA-based PCR assays.

Identification of citrus immune regulators involved in defence against Huanglongbing using a new functional screening system.

Huanglongbing (HLB) is the most devastating citrus disease in the world. Almost all commercial citrus varieties are susceptible to the causal bacterium, Candidatus Liberibacter asiaticus (CLas), which is transmitted by the Asian citrus psyllid (ACP). Currently, there are no effective management strategies to control HLB. HLB-tolerant traits have been reported in some citrus relatives and citrus hybrids, which offer a direct pathway for discovering natural defence regulators to combat HLB. Through comparative analysis of small RNA profiles and target gene expression between an HLB-tolerant citrus hybrid (Poncirus trifoliata × Citrus reticulata) and a susceptible citrus variety, we identified a panel of candidate defence regulators for HLB-tolerance. These regulators display similar expression patterns in another HLB-tolerant citrus relative, with a distinct genetic and geographic background, the Sydney hybrid (Microcitrus virgata). Because the functional validation of candidate regulators in tree crops is always challenging, we developed a novel rapid functional screening method, using a C. Liberibacter solanacearum (CLso)/potato psyllid/Nicotiana benthamiana interaction system to mimic the natural transmission and infection circuit of the HLB complex. When combined with efficient virus-induced gene silencing in N. benthamiana, this innovative and cost-effective screening method allows for rapid identification and functional characterization of regulators involved in plant immune responses against HLB, such as the positive regulator BRCA1-Associated Protein, and the negative regulator Vascular Associated Death Protein.

Codon Usage Bias Analysis of Citrus tristeza Virus: Higher Codon Adaptation to Citrus reticulata Host.

Citrus tristeza virus (CTV), a member of the aphid-transmitted closterovirus group, is the causal agent of the notorious tristeza disease in several citrus species worldwide. The codon usage patterns of viruses reflect the evolutionary changes for optimization of their survival and adaptation in their fitness to the external environment and the hosts. The codon usage adaptation of CTV to specific citrus hosts remains to be studied; thus, its role in CTV evolution is not clearly comprehended. Therefore, to better explain the host⁻virus interaction and evolutionary history of CTV, the codon usage patterns of the coat protein (CP) genes of 122 CTV isolates originating from three economically important citrus hosts (55 isolate from Citrus sinensis, 38 from C. reticulata, and 29 from C. aurantifolia) were studied using several codon usage indices and multivariate statistical methods. The present study shows that CTV displays low codon usage bias (CUB) and higher genomic stability. Neutrality plot and relative synonymous codon usage analyses revealed that the overall influence of natural selection was more profound than that of mutation pressure in shaping the CUB of CTV. The contribution of high-frequency codon analysis and codon adaptation index value show that CTV has host-specific codon usage patterns, resulting in higheradaptability of CTV isolates originating from C. reticulata (Cr-CTV), and low adaptability in the isolates originating from C. aurantifolia (Ca-CTV) and C. sinensis (Cs-CTV). The combination of codon analysis of CTV with citrus genealogy suggests that CTV evolved in C. reticulata or other Citrus progenitors. The outcome of the study enhances the understanding of the factors involved in viral adaptation, evolution, and fitness toward their hosts. This information will definitely help devise better management strategies of CTV.

Noemi Controls Production of Flavonoid Pigments and Fruit Acidity and Illustrates the Domestication Routes of Modern Citrus Varieties.

In citrus, the production of anthocyanin pigments requires the activity of the transcriptional activator Ruby. Consequently, loss-of-function mutations in Ruby result in an anthocyaninless phenotype [1]. Several citrus accessions, however, have lost the ability to produce these pigments despite the presence of wild-type Ruby alleles. These specific mutants have captivated the interest of botanists and breeders for centuries because the lack of anthocyanins in young leaves and flowers is also associated with a lack of proanthocyanidins in seeds and, most notably, with an extreme reduction in fruit acidity (involving about a three-unit change in pH). These mutants have been defined collectively as "acidless" [2-4]. We have identified Noemi, which encodes a basic helix-loop-helix (bHLH) transcription factor and which controls these apparently unrelated processes. In accessions of Citron, limetta, sweet lime, lemon, and sweet orange, acidless phenotypes are associated with large deletions or insertions of retrotransposons in the Noemi gene. In two accessions of limetta, a change in the core promoter region of Noemi is associated with reduced expression and increased pH of juice, indicating that Noemi is a major determinant of fruit acidity. The characterization of the Noemi locus in a number of varieties of Citron indicates that one specific mutation is ancient. The presence of this allele in Chinese fingered Citrons and in those used in the Sukkot Jewish ritual [5] illuminates the path of domestication of Citron, the first citrus species to be cultivated in the Mediterranean. This allele has been inherited in Citron-derived hybrids with long histories of cultivation.

Changes in Anthocyanin Production during Domestication of Citrus.

Mandarin (Citrus reticulata), citron (Citrus medica), and pummelo (Citrus maxima) are important species of the genus Citrus and parents of the interspecific hybrids that constitute the most familiar commercial varieties of Citrus: sweet orange, sour orange, clementine, lemon, lime, and grapefruit. Citron produces anthocyanins in its young leaves and flowers, as do species in genera closely related to Citrus, but mandarins do not, and pummelo varieties that produce anthocyanins have not been reported. We investigated the activity of the Ruby gene, which encodes a MYB transcription factor controlling anthocyanin biosynthesis, in different accessions of a range of Citrus species and in domesticated cultivars. A white mutant of lemon lacks functional alleles of Ruby, demonstrating that Ruby plays an essential role in anthocyanin production in Citrus Almost all the natural variation in pigmentation by anthocyanins in Citrus species can be explained by differences in activity of the Ruby gene, caused by point mutations and deletions and insertions of transposable elements. Comparison of the allelic constitution of Ruby in different species and cultivars also helps to clarify many of the taxonomic relationships in different species of Citrus, confirms the derivation of commercial varieties during domestication, elucidates the relationships within the subgenus Papeda, and allows a new genetic classification of mandarins.

Long-Term Field Evaluation Reveals Huanglongbing Resistance in Citrus Relatives.

Citrus huanglongbing (HLB) is a destructive disease with no known cure. To identify sources of HLB resistance in the subfamily Aurantioideae to which citrus belongs, we conducted a six-year field trial under natural disease challenge conditions in an HLB endemic region. The study included 65 Citrus accessions and 33 accessions belonging to 20 other closely related genera. For each accession, eight seedling trees were evaluated. Based on quantitative polymerase chain reaction analysis of the pathogen titers and disease symptoms, eight disease-response categories were identified. We report two immune, six resistant, and 14 tolerant accessions. Resistance and tolerance observed in different accessions may be attributed to a multitude of factors, including psyllid colonization ability, absence of pathogen multiplication, transient replication of the bacterium, lack of pathogen establishment in the plant, delayed infection, or recovery from infection. Most citrus cultivars were considered susceptible: 15 citrons, lemons, and limes retained leaves in spite of the disease status. Resistance and high levels of field tolerance were observed in many noncitrus genera. Disease resistance/tolerance was observed in Australian citrus relative genera Eremocitrus and Microcitrus, which are sexually compatible with citrus and may be useful in future breeding trials to impart HLB resistance to cultivated citrus.

Detection of citrus huanglongbing-associated 'Candidatus Liberibacter asiaticus' in citrus and Diaphorina citri in Pakistan, seasonal variability, and implications for disease management.

We report the detection of the huanglongbing (HLB)-associated bacterium 'Candidatus Liberibacter asiaticus' from both plants and insects in Pakistan and the seasonal variability in the numbers of 'Ca. L. asiaticus'-positive psyllid vector, Diaphorina citri. Our studies showed that 'Ca. L. asiaticus' was detectable from trees in areas with maximum temperatures reaching nearly 50°C (average maximum of 42°C). However, the bacterium was present at very low levels in psyllids both in summer (June to August) and autumn (September to November) in contrast to reports from Florida, where the bacterium was detectable at very high levels during October to November. We hypothesize that hot summer temperatures in Pakistan may interfere with acquisition and replication of 'Ca. L. asiaticus' in psyllids and may lead to dead or non transmissible 'Ca. L. asiaticus' in plants. Psyllid counts were very low in both summer and winter, showed a population peak ('Ca. L. asiaticus'-positive vectors) in spring, and showed a larger peak ('Ca. L. asiaticus'-free psyllids) in autumn. Natural thermotherapy during hot summers and a low vector population during environmental extremes may have played a major role in long-term survival of the citrus industry in Pakistan. These results may be useful in developing management strategies for U.S. citrus industries in Texas and California.

A six nuclear gene phylogeny of Citrus (Rutaceae) taking into account hybridization and lineage sorting.

BACKGROUND: Genus Citrus (Rutaceae) comprises many important cultivated species that generally hybridize easily. Phylogenetic study of a group showing extensive hybridization is challenging. Since the genus Citrus has diverged recently (4-12 Ma), incomplete lineage sorting of ancestral polymorphisms is also likely to cause discrepancies among genes in phylogenetic inferences. Incongruence of gene trees is observed and it is essential to unravel the processes that cause inconsistencies in order to understand the phylogenetic relationships among the species. METHODOLOGY AND PRINCIPAL FINDINGS: (1) We generated phylogenetic trees using haplotype sequences of six low copy nuclear genes. (2) Published simple sequence repeat data were re-analyzed to study population structure and the results were compared with the phylogenetic trees constructed using sequence data and coalescence simulations. (3) To distinguish between hybridization and incomplete lineage sorting, we developed and utilized a coalescence simulation approach. In other studies, species trees have been inferred despite the possibility of hybridization having occurred and used to generate null distributions of the effect of lineage sorting alone (by coalescent simulation). Since this is problematic, we instead generate these distributions directly from observed gene trees. Of the six trees generated, we used the most resolved three to detect hybrids. We found that 11 of 33 samples appear to be affected by historical hybridization. Analysis of the remaining three genes supported the conclusions from the hybrid detection test. CONCLUSIONS: We have identified or confirmed probable hybrid origins for several Citrus cultivars using three different approaches-gene phylogenies, population structure analysis and coalescence simulation. Hybridization and incomplete lineage sorting were identified primarily based on differences among gene phylogenies with reference to null expectations via coalescence simulations. We conclude that identifying hybridization as a frequent cause of incongruence among gene trees is critical to correctly infer the phylogeny among species of Citrus.

In vivo generated Citrus exocortis viroid progeny variants display a range of phenotypes with altered levels of replication, systemic accumulation and pathogenicity.

Citrus exocortis viroid (CEVd) exists as populations of heterogeneous variants in infected hosts. In vivo generated CEVd progeny variants (CEVd-PVs) populations from citrus protoplasts, seedlings and mature plants, following inoculation with transcripts of a single CEVd cDNA-clone (wild-type, WT), were studied. The CEVd-PVs population in protoplasts was heterogeneous and became progressively more homogeneous in seedlings and mature plants. The infectivity and pathogenicity of selected CEVd-PVs was evaluated in citrus and herbaceous experimental hosts. The CEVd-PVs U30C, G128A and U182C were not infectious; G50A and 108U+ were infectious but reverted back to WT and 62A+, U129A and U278A were infectious, genetically stable and more severe than WT. The 62A+ and U278A and U129A accumulated at higher levels than WT in protoplasts and seedlings respectively. The effect of specific mutations on the predicted secondary structure of the CEVd-PVs' RNA coupled with the infectivity and replication studies suggested complex structure-to-function relationships for CEVd.

Nucleotide sequence and genome organization of Dweet mottle virus and its relationship to members of the family Betaflexiviridae.

The nucleotide sequence of Dweet mottle virus (DMV) was determined and compared to sequences of members of the families Alphaflexiviridae and Betaflexiviridae. The DMV genome has 8,747 nucleotides (nt) excluding the 3' poly-(A) tail. DMV genomic RNA contains three putative open reading frames (ORFs) and untranslated regions of 73 nt at the 5' and 541 nt at 3' termini. ORF1 potentially encoding a 227.48-kDa polyprotein, which has methyltransferase, oxygenase, endopeptidase, helicase, and RNA-dependent RNA polymerase (RdRP) domains. ORF2 encodes a movement protein of 40.25 kDa, while ORF3 encodes a coat protein of 40.69 kDa. Protein database searches showed 98-99% matches of DMV ORFs with citrus leaf blotch virus (CLBV) sequences. Phylogenetic analysis based on the RdRP core domain revealed that DMV is closely related to CLBV as a member of the genus Citrivirus. DMV did not satisfy the molecular criteria for demarcation of an independent species within the genus Citrivirus, family Betaflexiviridae, and hence, DMV can be considered a CLBV isolate.